H = fir1(30, 0.1, 'low');
rtfilter([1],H);
i=1:100;
fs=16000;f1=200;f2=1000;f3=5000;
sin1=sin(2*pi*i*f1/fs);sin2=sin(2*pi*i*f2/fs);sin3=sin(2*pi*i*f3/fs);
sintot=(sin1+sin2+sin3)/3;
y=[];
sum=0;
xx=tic;
for j=1:1000;
for i=1:length(sintot)
    %tic;
    y(i)=rtfilter(sintot(i));
    %x=toc;
    sum=sum+x;
end
end
toc(xx);
sum=sum/100000